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Abstract 

We study some connections between the random moment problem 
and the random matrix theory. A uniform draw in a space of moments 
can be Ufted into the spectral probability measure of the pair (A, e) 
where A is a random matrix from a classical ensemble and e is a fixed 
unit vector. This random measure is a weighted sampling among the 
eigenvalues of A. We also study the large deviations properties of this 
random measure when the dimension of the matrix increases. The 
rate function for these large deviations involves the reversed Kullback 
information. 

Key words: Random matrices, unitary ensemble, Jacobi ensemble, spectral 
measure, canonical moments, large deviations. 
AMS classification: 15A52, 60G57, 60F10. 

1 Introduction 



The aim of this paper is to establish a strong link between random mo- 
ment problems ([5],[l5]) and random matrices. 

As a matter of fact, in the last decade, emphasis has been put on the 
asymptotic behavior of large random matrices. Such a study was first mo- 
tivated by problems arising in theoretical physics ([29], [H]). Generally, the 
distribution of these random matrices is characterized by some invariance 
properties. In the GOEjv (resp. GUEjv), the matrix is N x N symmetric 
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(resp. Hermitian) and, except for this constraint, has independent Gaussian 
entries. In the CUEat, A]\f has the Haar distribution on U(A^) (the group of 
N X N unitary matrices). We refer to [1] and f2T] for a general overview on 
the subject. In these examples, the distribution of both the random eigenval- 
ues Ai, . . . , Aat and eigenvectors are precisely known. These eigenvalues are 
(stochastically) independent of the eigenvectors and the matrix of the nor- 
malized eigenvectors is Haar distributed. The main object of interest is the 
empirical spectral distribution (ESD) of An, which is the discrete measure 

"'"'-wX^^v (1) 

When An is drawn randomly, the asymptotic properties of the sequence of 
random probability measures {h'^^^)N are obtained by a moment method 
using 

J a;V^)(dx) = itr(A^), (fc > 1) , (2) 

and {h'^^^)N converges towards a deterministic probability measure. Namely, 
for the GOE (resp. GUE) the limit is the semicircle distribution while for 
the CUE, it is the uniform measure on the unit circle T. To give the rate of 
convergence of the sequence {u^^'')n, it is interesting to know if it satisfies 
the large deviation principle (LDP). This problem has been first studied for 
the GOE in the pioneer work of Ben Arous and Guionnet (|3j). It was further 
investigated for the CUE (see [20]) and many other distributions (see |13j). 
In these LDPs the speed is N'^ and the rate function may be obtained from 
the Voiculescu entropy ([21J). 

In this paper, we consider another approach to unitary matrices. Instead 
of considering the classical empirical measure u^^^ defined in ([1]), we study 
the random spectral measures (see (jl]) below). Every unitary operator A on a 
Hilbert space with a unit cyclic vector e is (by the spectral theorem) unitarily 
equivalent to the multiplication by z on L^(T, /x) for some probability measure 
/i, which is uniquely determined by its moments. If we endow U(A^) with the 
normalized Haar measure Au(Ar), the first vector of the standard basis ei is 
a.s. cyclic for An- We thus define the random spectral measure /xi^^ on T 
related to the pair {An, Ci). It is defined by its moments: 

[ ;^Vf)(dz) = (ei,AW, (^>1)- (3) 
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The ESD u^^^ and the spectral measure yU„ are widely different. This 
can be seen in terms of moments. In both cases, we consider successive 
powers of A^, but the moments of fii^^ are the first entries given by ([3]) and 
those of u^^^ are the normalized traces given by ([2]). Alternatively, a simple 
diagonalization of An gives: 

/^f^-EVe-.. (4) 

where the w^-'s are the square moduli of the top entries of the normalized unit 
eigenvectors and exp(z6'j) = Xj are the (a.s. different) eigenvalues of An- It 
is now clear that this measure carries more information on An than z/^^\ 
Ideed, it involves weights depending on the eigenvectors of this matrix. In 
this paper, we are concerned with a precise study of the random probability 

(N) 

measure /i„ . 

For a probability measure ^ on T, two encodings are interesting. First, 
its Verblunsky coefficients appear in the Schur recursion formula satisfied 
by the sequence of orthogonal polynomials in L^(^) (see for example jM] . 
|35j). Secondly, its canonical moments are involved in statistical applications. 
They are inductively defined : the (A; + l)-th canonical moment is the relative 
position of the {k + l)-th moment consistent with the k first ones (see [TT] 
Chapter 9 for an overview on canonical moments). 

Two remarkable papers motivated our work. On the one hand, Killip and 
Nenciu |26] proved that in the CUEat model, the Verblunsky coefficients of 

^ are independent and have known distributions (namely complex beta 
distribution). On the other hand, Lozada [27| has studied the distribution 
of the N first canonical moments for special random measures, namely the 
measures having their first moments uniformly distributed. Surprisingly, 
the distributions found by Killip and Nenciu and by Lozada are the same. 
We can explain this last result by the fact that Verblunsky coefficients and 
canonical moments coincide (as pointed out by several authors, see Simon 
[53] p. 439). Using this observation we show in Theorem 14 . 1 1 that the vector 
of moments of /il^'' under CUEtv is also uniformly distributed, giving a (new) 
enlightening connection between random moment theory and random matrix 
theory. 

Extending this method, we give the precise distributions of moments 
and canonical moments in other circular ensembles and in the so-called (3- 
ensembles (or log-gases). Furthermore, we consider real matrices. Beginning 
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with the special orthogonal group of order 2N SO(2A^), the symmetry prop- 
erty of spectral measures allows to consider their projections on [—1, 1], which 
leads naturally to the family of Jacobi ensembles. We prove a lift of a uni- 
form random point in the space of moments into the spectral measure of a 
random element of the ensemble SO(2iV)/U(iV) (Theorem I4.10p . 

The common feature of all these models is that the weights and the sup- 
port points of jjL^^^ are independent. It is a nice framework to study the 
LDP, with two main differences with the above mentioned LDP for ESD : 
the (random) weighting both slows down the speed from N"^ to and (con- 
sequently) changes the rate function to the reversed Kullback information. 
This representation allows in particular to recover results proved by using 
the independence of canonical moments in [12] and [2Z!- As a consequence, 
we see that yul^'' converges weakly in probability to the same deterministic 
limits as the ESD z/(^). 

Some results for ensembles with non compactly supported spectral mea- 
sures can be found in [2] for sample covariance matrices, and in [16] for the 
GUE. 

The paper is organized as follows. In Section [2], we give some notations 
and list some distributions of frequent use. In Section [31 we explain the set- 
up of moments and canonical moments in the complex and real case. Section 
m is devoted to the distribution of random spectral measures in different 
models. Finally, in Section [5], we state LDPs for the families of random 
spectral measures. 

Notice that a few months after dropping the first version of the present 
paper on Arxiv, we have been aware of the work of Birke and Dette |1] 
on the asymptotic behavior of the roots of random orthogonal polynomials 
associated to random moments. Their main result is the computation of the 
root distribution when the moments are uniformly distributed. It appears to 
be a particular case of our more general Lemma 14.71 and Theorem 14. 8[ 

2 Notations and some useful distributions 

Let D be the open unit disk in C and T = 5©, the unit circle parameter- 
ized hy z = e*^, with 9 G [— 7r,7r). Let /i be a probability measure on T. ix 
is said to be trivial if it supported on a finite set, and nontrivial otherwise. 
Further, the Hermitian product on L^(T, /i) is defined by (/, g) := Jjfg d/i. 
Let (G, Q) be any measurable space, we denote by A4i{G) (resp. M.{G)) the 
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set of all probability measures (resp. positive measures) on G. If G Ai\{G) 
and / is integrable, we write sometimes for / d/x. We recall now some 
special useful distributions that are used later. For A; > 1, we set 

Sk {{xi, • • • , Xfe) : > 0, (i = 1, • • • , A;), xi H h Xfe = 1} 

:= {(xi, • • • , Xfc) : > 0, (i = 1, • • • ,k), Xl^ h < 1} . 

Obviously, the mapping (xi, • • • , Xk+i) i— > (xi, • • • , Xk) is one to one from the 
simplex 5^+1 onto . 

For > 0, J = 1, . . . , /c + 1, the Dirichlet distribution Dir(ai, • • • , a^+i) 
on Sk+i has the density 

r(ai H hQfc+l) afe+i-l .^x 

r(aO---r(a,+o 

with respect to the Lebesgue measure on Sk+i- When ai = • • • = a^+i = a > 
0, we denote the Dirichlet distribution by Diik+i^o). If a = 1 we recover the 
uniform distribution on Sk+i- 

Pushing the Dirichlet distribution under the previous mapping, we get a 
probability measure Dir(ai, . . . , a^; ak+\) on with density : 

r(ai H h gfc+l) _ _ _ , Na,+i-l 

r(aO---r(a,+o '^'^ "^^^ 

If A; = 1, Dir(a; h) has the Beta(a, h) density on [0, 1] : 

r(a + 6) , 

r(a)r(6) ^ ^ ■ 

The very particular case Beta(l/2, 1/2) is the so-called arcsine law. Some- 
times we need the distribution obtained by pushing forward Beta(a, h) under 
the mapping x ^-^ 2x — 1. It is the distribution Beta5(6, a) on (—1, 1) having 
density 

r(a)r(&) + ■ 

We use also a complex version of the beta distribution. For r > — 1 let 77,. be 
the probability density on C defined by 

nr{z) ^ (1 - \A'Y ' e D). (6) 



It is obviously the density oi X = e^^ \/B where U is uniform on [0, 2tc] and 
B is Beta(l, r) distributed. 

To end this section, let us recall the classical relation between the Dirichlet 
and Gamma distributions. A Gamma variable ja with parameter a > has 
density 

r(a) 

It is well known that if i = 1, ■ ■ ■ , r are independent and if yi = ■jbi {h > 0) 
then 

— \ : : ] =BiT{bi,--- ,br) (7) 

^yi + --- + yr yi + --- + yrj 

and this variable is independent of yi + ■ ■ ■ + i/r- 



3 Verblunsky coefficients, canonical moments 
and uniform probability 

We now present the notion of canonical moments of a measure. We 
successively focus on the complex case (unitary circle) and on the real case 
(compact interval). Finally we recall that the uniform probability on moment 
spaces corresponds to a nice distribution on the space of canonical moments. 



3.1 The complex case : Verblunsky coefficients 

All the material of this subsection comes from [31] Section 1 or [53] Sec- 
tions 2 and 3. We recall here the connection between moments of a prob- 
ability measure on T and Verblunsky coefficients built through orthogonal 
polynomials. 

Let fi be an arbitrary nontrivial probability measure on T. The functions 
\^ z ^ z 5 ' ' ' are linearly independent in L^(T, yu). Following the Gram-Schmidt 
procedure we define the infinite sequence ($j)j>i of monic orthogonal poly- 
nomials. More precisely, ^o{z) = 1 and ^j{z) is the projection of z^ onto 
{l,---,^^-i}^,forj>l. 

If /i G M.i{T) has finite support {^i, ■ ■ ■ ^zk} {K different points), we 
still define $j in the same way for j = 1, . . . , — 1. Besides, we define 
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as the unique monic polynomial of degree K such that = i.e. 

K 

^k{z) = \{{z - z,) . (8) 

By convention, in the nontrivial case we set = oo, and then the sentence 
"for every j < K" (or "for every j < K + 1") mean "for every finite j". 

Some useful polynomials associated to the sequence {^j)j<:K are the re- 
versed (or reciprocal) polynomials. They are defined by $o(^) = -'- 

^*^iz) = z^^JiT/^ . (9) 

In the next proposition, we define the Verblunsky coefficients. They are 
sometimes called Schur, Szego or Geronimus coefficients or even reflection 
coefficients. Surprisingly, these coefficients also appears as central quantities 
in the theory of moment problems (see for example [H]), where they are 
called canonical moments. 

Proposition 3.1 (Szego) For j < K +1, we define the Verblunsky coeffi- 
cient of order j by setting cj := — $j(0). The sequence {^j)j^K+i satisfies 
the recursion 

$0 = 1 , = ^$,-1 - c,<l>^*„i il<j<K+l). 

Notice that if fi is nontrivial, cj G © for every j > 0. If K < oo, then 
Cj e D for 1 < J < /sT - 1 and ck e T (see |3l] Theorem 1.5.2). In the 
sequel, when we need to stress the dependence of the Verblunsky coefficients 
on the underlying measure /i, we write Cj{fi). A theorem due to Verblunsky 
claims that the correspondence between a probability measure on T and the 
sequence of its coefficients is one to one ([33] Theorem 1.7.11). For > 1, 
set 

Mj, = \( [ z'fiidz)] : /i e Mi{T) \ . (10) 

I V^T / l<j<N J 



Proposition 3.2 ([34] pp. 60 and 218, [TT] p.269) Let{ti,--- ,tN) e int M' 
the range of the [N + l)th moment 

tN+i = [ z'^^'r^idz) 
Jt 
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as 7] varies over all probability measures having (ti, ■ ■ ■ ,tj\f) as their N first 
moments, is a disk centered at some point s^r (depending onti,--- ^t^) with 
radius 

N 

= n (1 - i^^f ) ^ 

The relative position is 
and it is exactly cn+i ■ 

This striking connection of two notions (related position versus recursion 
coefficient) has been proved in two opposite ways: starting from the recursion 
to get relative position [34J p. 60), or alternatively starting from the moments 
and their relative positions and computing the recursion [33] p. 218, 
p. 269. For an historical account, see [33] P-10 and p. 221. 



3.2 Canonical moments : Real case 

Following the last interpretation of canonical moments as relative posi- 
tions of moments in M^, it is possible to build similar quantities when the 
set of integration T is replaced by a compact interval [a, b] {a < b). Let us 
briefly recall the exact definition of canonical moments in this frame. We 
refer to the excellent book of Dette and Studden [TT] for a complete overview 
on the subject. To begin, for > 1, we denote by m]J''''' the A^-th moment 
space generated by probability measures on [a, b\. 



x''fx{dx),k = l,...,N^ : ^eMi{[a,b])^ 



Given m-' := (mi,-- - ,mj) E int Mj , (j > 1), we first define the 
extreme values, 



m++i(m^) 



= max |m G M : (mi, - - - , ruj, m) G iVfj^i | (11) 
= min |m G M : (mi, - - - ,mj,m) G Mj'^J^ | . (12) 
For k > i > 1, the i-th canonical moment is defined recursively as 



m^.+i(m^) 



c.(m^)^c.(m-):= '"• " (13) 
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A quite nice property of canonical moments is that they are invariant on 
any affine one to one mapping transforming the support of the underlying 
measures (see for example plj). So that, we may restrict ourselves to the 
special case of M^^'^K 

3.3 Uniform probability on moment spaces 

In this subsection, we endow the moment sets and Mn'^^ (for fixed 
n) with the uniform probability and recall that in these cases the canonical 
moments previously defined have very interesting properties (see Lemmas 13 .41 
and 13.31 below). 

As pointed in Section 13.11 for fixed there is a one to one mapping 
between moments and canonicals moments. Indeed, we may define a one to 
one mapping kJ^ 

kI : int Mjf 

(ti,...,t7v) ^ {Ci,...,Cn). (14) 

In the real case, the relation f|T3l) defines a one to one mapping k^^'^^ 

4^' : int MjJ'^l (0, 1)^ 

(tl,...,tAr) ^ (Ci, . . . ,C7v) . (15) 

which is triangular and bijective. 

The two following lemmas give the canonical moment distribution when 
the moments are uniformly drawn. 

Lemma 3.3 (Lozada) Endowing int Mjj with the uniform distribution is 
equivalent to the N first canonical moments (ci, . . . ,cn) being independent 
in such a way that Cj has density rjiq_j. 

Lemma 3.4 (Chang, Kemperman, Studden) Endowing int m]^'^^ with 
the uniform distribution is equivalent to the N canonical moments (ci, . . . , cat) 
being independent in such a way that Cj is Beta(A^ — j + 1, A^ — j + 1) dis- 
tributed. 

These two lemmas have been the starting point for the investigation on the 
asymptotic behavior of the randomized sets of moments m]^'^^ and Mjf (see 
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[5], [15], [2Z] and [ID])- Here, these lemmas are useful to obtain some very 
nice elementary properties of random measures built on eigenvalues of some 
classical random matrix models. We develop these results in the next sec- 
tions. 

4 Random spectral measures and their mo- 
ments 

In this Section, we first define the spectral measure associated with the 
pair {A, e) where A is a unitary matrix and e a cyclic vector. Then we 
randomly draw a matrix and study the distribution of the associated random 
spectral measure and of its canonical moments. We focus on various classical 
distributions popular in the random matrix paradigm. Moreover we extend 
the class of such random measures leading to log-gases models, which in turn 
may be lifted into matrix ensembles. In particular, we emphasize the different 
ways to get uniform distribution on the space of moments. Our main results 
in this section are Theorem 14. Theorem 14.81 and Theorem 14.101 



4.1 Spectral measures associated with a unitary ma- 



Let us consider with its canonical basis (ci, . . . ,ej\f) and 1]{N) the 
group of N X N unitary matrices. For A G U(A^), let Aj, i = 1, . . . , N be the 
eigenvalues of A. We may write 



where D := diag(Aj)j=i,.,.,Ar is diagonal and 11 := {TTij)ij=i,...,N is unitary. 
Let us assume that ei is cychc, i.e. that span(ei, Aei, . . . , A^'^ei) has rank 
A^. Following [3S], we consider the spectral measure /xi^^ associated with 
the pair {A,ei) i.e. having, for any G IM, moment of order k equal to 
(ei,y4'^ei). Obviously, /il^'' is trivial and supported by the eigenvalues of A. 
More precisely. 



it 

where := IvTifcp, k = 1, . . . , N are the square moduli of the top entries of 
the normalized eigenvectors. 



trix 



A = UDU*, 



(16) 




(17) 
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4.2 Circular ensembles 



Here are recalled three families of classical ensembles of random matrices: 

• CUE^ (circular unitary ensemble) : U(iV) equipped with its normalized 
Haar measure \u(n) 

• COEtv (circular orthogonal ensemble) : set of symmetric unitary NxN 
matrices. Every element 5* can be written S = g^g with g G U(A^) 
hence COEat = U(A^)/0(A^). The probability measure is the pushfor- 
ward of Xi](n) under the projection mapping U(A^) U(A^)/0(A^). 

• CSEjv (circular symplectic ensemble) : set of 2N x 2N self-dual unitary 
matrices. Recall that the dual of a matrix H is := JH^J^ with 




and that a 2N x 2N matrix k is symplectic if it satisfies kJk = J. 
Every element 5* can be written S = g^g with g G U(2A^), hence 
CSEjv = U(2A^)/USp(2A^), where USp(2iV), the set of unitary sym- 
plectic matrices, is the invariant set of the involution g \—>- {g^)~^. The 
probability measure is the pushforward of Au(2Ar) under the projection 
mapping IS{2N) U(2iV)/ USp(2A^). 

All these ensembles are considered in the two following subsections. 



4.2.1 The circular unitary ensemble 

Let us first notice that with probability 1, the vector ei is cyclic. If 
A G U(A^), all its eigenvalues have unit modulus and in the decomposition 
(fT6|) we may write D = diag (e*^^, ■ ■ ■ ,e^^'^Y Further, the matrix 11 in (fT6!) 
is not uniquely determined. But, it is well known that it can be chosen such 
that, if A is Haar distributed then H is also Haar distributed and independent 
of D. In this case, the joint law of eigenvalues arguments (on [0, 27r]^) is 

CUE^(d^i, ■ ■ ■ ,d0N) := |A(e^^S ■ ■ ■ ,e^^-)|2 d^i ■ ■ -d^^ (18) 

where A is the Vandermonde determinant 

A(xi, ■ ■ ■ ,XAr) = Yl (Xj-Xk). (19) 

^<j<k<N 
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Furthermore, since 11 is Haar distributed, its first row is uniformly distributed 
on the Euchdean sphere of C^. So that, setting := |vri ^p, {k = 1, . . . , N) 
the vector (wi, • ■ ■ is uniformly distributed on Sn (see for example 

Proposition 3.1 of [26]) i.e. 

(wi,...,w;v) = Dir;v(l) (20) 

(see ([5]) for the definition). 

So, from f[T7|l we have built the random probability measure on T: 

N 

-"l'^^ = E^'^^e-., (21) 
k=l 

supported by the eigenvalues of the matrix A with an independent system of 
weights. 

Up to our knowledge, the following result was not known previously. 
Theorem 4.1 Under CUEtv, the distribution of ( Lz^ai'^\dz)) is 

\ / l<j<N-l 

uniform on Mj^^^. 

Proof: Let {ci{fii^^), ■ ■ ■ , CAr(/il^^)) be the vector of first canonical mo- 
ments of fii^\ In Proposition 3.3 of [26] (see also [36] Section 11) it is proved 
that the components of this vector are independent. Furthermore, 

CUE^(c,(^f )) G d^) = Vn-,-i{z) dz, j = 1,--- ,N-1, (22) 

and ctv is uniform on T (since the support of finite set, the last 

canonical moment belongs to T). Considering now the pushforward of these 
distributions under (/t]^_i)~^ the measurable one to one mapping from D^~^ 
to int Mjf_i which map canonical moments to moments we may conclude 
using Lemma [3^31 ■ 



For the special unitary group SU(A^) we also have. 

Corollary 4.2 Under the normalized Haar measure Xsv(n) on E>IJ{N), the 
distribution of ( Jj z^ jj,i'^\dz)j is uniform on Mjj_-^^. 

V / 1<J<A''— 1 
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Proof: The joint eigenvalue distribution on ^ of the normahzed Haar 
measure on §U(A^) is 

^ - |A(e^^S--- ,e*^-)p d^^i ■ ■ ■ d^^^._i 



.N- 



with 9iy = —{9i+- ■ ■+6'7v-i)(niod27r) (see Notice that it is the density of 
(6*1, . . . , Om-i) under CUE at conditioned on 9n = —{Oi + ■ ■ ■ + 6'7v_i)(mod27r). 
Moreover we have exp {i{9i + ■■■ + On)) = (-l)^$jv(0) = (-l)^+^cjv, so 
that the conditioning set is {cat = (—1)^^^}. We have seen in the proof 
of Theorem 14.11 that Ci, • • ■ , cat are independent under CUEtv- Therefore, 
Ci, ■ ■ ■ , CAr_i remain independent when conditioning on {cjy = (— 1)^"*"^}. 
Moreover, their distributions are not affected by this conditioning. Using 
Lemma [3.31 we deduce that the conclusion of the Corollary holds true. ■ 



4.2.2 The circular /9-models 

For the COE and CSE we have also independence between the weights 
and the support points. Moreover, the distributions are given by formulas 
analogous to f|T8l) and fl20l) . Actually it is convenient (and now classical) to 
define a class of distributions with a continuous parameter f3 > 0, extending 
the cases (3 = 1 (COE), f3 = 2 (CUE) and /3 = 4 (CSE). 

Definition 4.3 For P > 0, let C/5Ejv the distribution on x Sn such that 

1 ) The joint law of {6i, . . . , Oj^j) is : 

C/3E^(d^i,--- ,d^^) :=C^(iV)-i|A(e^^S--- ,e^^")p d^i---d^^, 
where 

r (1 + 

(r(i + f)r 

2) The joint distribution of (wi, ■ ■ ■ , Wjv) on Sn is DirAr(/5/2). 

3) The variables 6i, - ■ ■ ,0n and wi, ■ ■ ■ , wat are independent. 
Setting 

N 

/xf) = 5^w,5,.., (23) 

k=l 



13 



PB] Killip and Nenciu proved in their Proposition 4.2 that under C/3E^ : 

1) the Verblunsky coefficients Cj = Cjd^i^^), {j = 1, . . . ,N) are indepen- 
dent, 

2) forj = l,...,iV-l, 

C/?E^(c,(/xf )) G dz) = dz , (24) 

3) the distribution of CAr(/il^'') is uniform on T. 

Actually, they built also an explicit random matrix model (namely a 
five-diagonal matrix) whose spectral elements (e*^^, . . . , e*^'^) and (wi, . . . , wat) 
satisfy the previous properties. 

4.3 SO{2N) and Jacobi ensembles 

In this subsection, we consider real random matrices. We ffist consider 
the subgroup SO(2A^) of U(2A^), for which spectral measures (on T) are 
symmetric (i.e. invariant by complex conjugation). This entails that the 
Verblunsky coefficients are real. We recall the results of [26] on the distribu- 
tions of these coefficients, for SO(2iV) equipped with the Haar measure and 
also for a three-parameter family of log-gases. Then, we project such spectral 
measures on R and, by tuning the parameters we obtain random probability 
measures on [0, 1], whose moments are uniform. 

4.3.1 §0(2iV) and extension 

If we provide §0(2A^) with the normalized Haar measure A§o{27V); then 
±1 are a.s. not eigenvalues and the spectral measure defined in f[T7l) is a.s. 
supported by pairwise conjugate complex numbers. It may be written as: 

N , 

/^r^ = Ey(^e-.+5e-.)- (25) 
k=l 

In that case, the distribution of (6*1, •• ■ , 6]^) has a density proportional to 

|A(2cos^i,--- ,2cos^jv)P (26) 

and the array of weights (w'^,-- - ,w'jv) is independent of (6'i, ■ ■ ■ ,0n) and 
satisfies 

(w'l, ■ ■ ■ ,w'^) = DirAr(l) , 
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(see [2Sj Proposition 3.4). 

A natural generalization is obtained replacing the square in fl26l) by an ex- 
ponent (3 as in Section [4.2.21 and adding a particular external potential: 

Definition 4.4 For a, b, (3 > let J{(3, a, b, N) be the distribution of{6i,--- , 9^) 

having a density proportional to 

N 

I A(cos ^1, ■ ■ ■ , cos On) f - cos (1 + cos (27) 

k=l 

In this generalization, the random canonical moments have some very re- 
markable properties (like for the C/5E^ model) as soon as the weights are 
suitably distributed: 

Proposition 4.5 (Killip-Nenciu J2B^ Prop. 5.3) 

1. For SO(2A^) equipped with the Haar measure \eiO){2N), the Verblunsky 
coefficients Ci{fi^^^), ■ ■ ■ , C2n-i are independent and satisfy: 



Moreover C2Ar(/ul^^'') = — 1. 

2. More generally, under the distribution J{(3,a,b, N) (g) DirAr(/?/2) on 
[0,27r]^ X Sn, the Verblunsky coefficients ci{iJ,^^^), ■ ■ ■ ,C2N-i{fJ'^^^) 
are independent and satisfy: 

Ck{^^^^>) = Beta, i f3 + a, f3 + b\ k odd 

/ (2N), W T. f2N-k-2^ , 2N-k\ , 

Cfc(/i; ') = Beta, I p + a + b, — - — pi k even, 



for k < 2N — 1. Moreover C2N{f^^^^) = ~1 



It is then natural to project the symmetric measure ni^^ on M. It is the 
motivation of the next subsection. 



15 



4.3.2 Jacobi ensembles 

Let fi he a. symmetric probability measure on T (i.e. invariant under 
complex conjugation) and let 7 = R{fi) be the pushforward of ^ by the 
mapping e*^ 1— > kt^sA^ go that for g continuous on [0, 1] 

I' 9ix) d7(x) = ^ 9 (^^^f^) Kdz) . (28) 
Applying that to a probability measure with finite support 

N , 
fe=l 

we get a probability measure 

7f)=^(^l^^) = E-'A. where X, = i±|^ , fc = l,...,iV. (29) 

k=l 

The correspondence between the Verblunsky /canonical coefficients of /j, and 
the canonical coefficients of 7 = R{fi) (in the sense of Subsection 13.21) is 

Cfc(7) = ^ + cM) (or. equivalently Cfc(/i) = 20^(7) - 1) . (30) 

This claim is proved by a geometrical argument in [32] (Theorem 13.3.1) and 
using Tchebychev polynomials in [TI] (Section 9.3.8). 

Definition 4.6 For a,b,P > 0, let J{P, a,b, N) the distribution on [0,1]^ 
with density proportional to 

N 

\A{xu...,XN)fl[xl\l-Xkr-\ 

k=l 

If we equip [0, 27r]^ with J(/5, a, b, N) then the joint distribution of Xi, ■ ■ ■ x„ 
is J(/9, a, b, N). We are now ready to claim : 

Lemma 4.7 Under the distribution J{/3, |, |, A^)(g)Dir7v(/3/2), the canonical 
moments 01(7^^^), ■ ■ ■ , C2Ar-i(7w^^) are independent and 

c,(7fV^^Beta(^^/?,^^/5) (A; < 2iV - 1) . (31) 
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Proof: From fl30|) and fl29|) we have Ck{'yi^^) = |[1 + CA;(/il^'')]. Applying 
Proposition 14.51 2) . we see that the new coefficients inherit independence and 
have the mentioned distributions. ■ 



Theorem 4.8 



1. If-fl, = Y.k=l^k^^k ^^^^ 



(xi, . . . , x;v, wi, . . . , w^) J(4, 1, 1, N) ® Dir^(2) , (32) 
then the distribution of (mi(7w^''), ■ ■ ■ , m2Ar-i(7i^'*)) is uniform on 



(xi, . . . , x^, Wo, . . . , w;v) = J(4, 1, 3, N) ® Dir(l, 2, ■ ■ ■ , 2, 2) , (33) 
then the distribution of (mi{'yi^^), ■ ■ ■ , m2N{li^^)) is uniform on M^^^ . 
3. If'yi^^ = wo(5o + Ylk=i ^kSxk + wat^i with 

(xi,...,x;v-i,wo,...,w^) J(4,3,3,iV-l)®Dir(l,2,-- - ,2,1),(34) 



then the distribution of [mi{'-fi^^) , ■ ■ ■ ,m2N-i{'yi^^)) is uniform on 

4. If'ji^^ = Ylk=i^kSxk +^N+iSi with 

(xu . . . , x;v, wi, . . . , wjv+i) = J(4, 3, 1, AT) ® Dir(2, 2, ■ ■ ■ , 2, 1) , (35) 

then the distribution of (mi{'yi^^), ■ ■ ■ , m2Ar(7w^'')) is uniform on M^j^^ 

To prove the theorem we need the following proposition on principal rep- 
resentations of finite moment sequences (see for example [23], [25] or [11] 
Defininition 1.2.10). These representations are related to the extreme values 
introduced in ffTTl) and ffT2D. 



Proposition 4.9 
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1. Let (mi, ■ ■ ■ , m2N^i) e Mf^\^ 

• Lower principal representation: there exists a unique system < 
Xi < . . . < xat < 1 and wi, . . . , wjv G (0, 1) with Yl!k=i '^k = ^ such 
that the probability measure z/_ := Yl^^i^k^Xk satisfies 



rrij = / x-'h'^(dx), j = 1 . . . 2N — 1. 
Jo 



Upper principal representation: there exists a unique system < 
Xi < ... 
such tha 
satisfies 



Xi < . . . < Xn-1 < 1 and wq, . . . , Wjv € (0, 1) with J2k=o = ^ 
such that the probability measure := wq^o + '^k=i ^k^x^ + '^N^i 



rrij = / x-'V+(dx), j = 1 . . . 2N. 
Jo 



2. Let (mi, ■ ■ ■ , mg^v) e int (M^^n^) 

• Lower principal representation: there exists a unique system < 
xi < . . . < xat < 1 and wq, . . . , wat G (0, 1) with XlfcLo ^fc = 1 such 
that the probability measure U- := wq^o + J2k=i "^k^x^ satisfies 

rrij = / xV_(dx), j = 1 . . . 2N. 
Jo 

• Upper principal representation: there exists a unique system < 
Xi < . . . < Xat < 1 anc? Wi, . . . , Wtv+i G (0,1) with J2k=i '^k = 
1 such that the probability measure := Ylk=i^kSxk + "^N+i^i 
satisfies 

rrij = / x-'u+{dx), j = 1 . . . 2N. 



Proof: [of Theorem US] 

• To prove the first point, take /3 = 4 in Lemma [4.71 and n = 2N — 1 in 
Lemma 13. 4[ 

• To prove the remaining parts of the theorem we use Proposition 14.91 
We only treat the point 2) since the proof of the other ones may be 
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tackled identically. From Proposition 14.91 (even case, lower representa- 
tion) the mapping between {xi, ■ ■ • , Xn, Wi, ■ ■ ■ , w^v) and (mi, ■ ■ ■ , m2N) 
is a diffeomorphism. Thus the probability density of the moments are 
completely known once we compute the Jacobian. Hence, we have to 
compute 



J = 


Xl 

xl 

xl 


Wl 

2wiXi 
Swi^i 


X2 
^2 
•^2 
^3 
•^2 


W2 
2W2X2 
3^2X2 


■ Xn 


Wat 

2\}]s[Xn 




•^1 


2Arwix2^-^ 


•^2 


2NM2xf'-^ ■ 


X2N 


2iVw^x2j^-^ 



An easy calculation yields 



J 




fc=i / 



Up to a normalizing constant, we recognize the density of J(4, 1, 3, A^)c 
,2,2). I 



Dir(l,2, 



From Lemma [331 the distribution of the canonical moments in the four cases 
of the previous theorem is completely known. In their Section 2, Killip and 
Nenciu [26j give a model of tridiagonal real matrices admitting these spectral 
characteristics: their spectral random measures have the same distribution 
as 7„. Here, we present for the first and fourth cases of Theorem 14.81 another 
interesting model issued from symmetric spaces. 

Theorem 4.10 Let U be a 2n x 2n random matrix such that U = g^g where 
g is Haar distributed on SO(2n). The spectral measure (i^^ ofU is symmetric 
and the distribution of (7/11(7^"''), ■ ■ ■ ,mn-i{'yi"'^)) is uniform on M^^^. 

Proof: Assume first that n = 2N. From Theorem 14.81 1. it is enough to 
check that this model induces the distribution J(4, 1, 1, A^) (g) DirAr(2). We 
follow the notation of |12j . 

Let $ be the mapping g G E>0(2n) 1— > g^g and let S(n) be its image. 
Since $~^(1) := K{n) = SO(2?2) fl USp(2n) we can see $ as a one to one 
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mapping from E>0{2n)/K{n) onto S(n). Actually K{n) is isometric to l]{n) 
via 

Let 5* be a random element of S(?7,) whose distribution is the pushforward of 
A§o(2n) under $. We have 

S = i{g)a{i{g))-^ 
where g and a are independent, g is Au(n) distributed and if n = 2N 

where A^ = diag(e*^% . . . , e*^'^) and the 9k s have a joint density proportional 

to 

N 

A(cos6'i, . . . , cos6'iv)'^ Y\ I sin6'fc| . 

k=l 

(see [12], proof of Theorem 2 (formulae 47 and 48)). In terms of Xk = 
(1 + cos6k)/2 this yields exactly the J(4, 1, 1, A^) distribution. Moreover, a 
simple computation gives 

TV 

(ei,S'^ei) = ^w^cos(j6'r) 

r=l 

with Wr = (ki.rP + \gi,N+r\'^) a^d (wi, . . . , w^v) is DirAr(2) distributed. 

If n = 2A^ + 1, similar considerations may be made to interpret formula 
( |35|) as the random spectral measure of §0(2n)/A'(n) when n = 2N + 1 (in 
that case 1 is double eigenvalue of S"). ■ 

5 Large deviations for random spectral mea- 
sures 

We present here the LDP for sequences of random measures defined in 
the previous section. The main result is Theorem 15.31 which is obtained as 
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a consequence of our generic Proposition 15.21 concerning weighted random 
measures. This proposition states that the LDP for empirical measures with 
speed A^^ can be used to derive the LDP with speed for measures with 
Dirichlet distributed weights. 

For the sake of completeness we briefly recall the LDP definition. Let {un) 
be a decreasing positive sequence of real numbers with lim^^oo = 0. 

Definition 5.1 We say that a sequence (Rn) of probability measures on a 
measurable Hausdorff space {G, B{G)) satisfies the LDP with rate function I 
and speed (uj/) if: 

i) I is lower semicontinuous (Isc), with values in IR"*" U {+00} . 

ii) For any measurable set A of G: 

— /(int A) < liminf Mat log-RAr(A) < lim sup mat log -Rat (A) < — /(clo A), 

where I{A) = inf^g^/(^) and int A (resp. clo A) is the interior (resp. 
the closure) of A. 

We say that the rate function I is good if its level set {x G G : I{x) < a} 
is compact for any a > 0. More generally, a sequence of G -valued random 
variables is said to satisfy the LDP if the sequence of their distributions 
satisfies the LDP. 

Hereafter, the space G is the set of all non negative measures (or probability 
measures) supported by a given compact set. We endow this set with the 
topology of the weak convergence. The rate function obtained in this paper 
is be the so-called reversed KuUback information. Let us recall that if P and 
Q be probability measures on G, the KuUback information between P and 
Q is defined by 

. , X f /log^dP ifP<Q and log^ eLi(P), 
JC{P\Q) = Ug dQ ^ ^'iQ ^ ^' (36) 

\^ +00 otherwise. 

For fixed Q, the nonnegative convex function }C{-\Q) is the rate function of 
the LDP for the empirical distribution of a sample according to Q (Sanov's 
theorem [9] p. 263). For fixed P, the function /C(P|-) is nonnegative convex, 
it vanishes only for Q = P, and we call it reversed Kullback information with 
respect to P. 
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5.1 LDP for random measures with Dirichlet weights 

Proposition 5.2 Let p > 1 and {C,k,N, ^ ^ k < A^}Ar>i be a triangular array 
of random variables belonging to a compact interval [xi,X2] and let 

1 ^ 

k=l 

Assume that the sequence (Cn) satisfies the LDP in A^i([xi)X2]) with speed 
{NP) with a good rate function I^. Assume further that has a unique 
minimum at v whose support is [xi-,X'2\- Let {\}k,N, 1 < ^ < another 
triangular array, independent of the first one, such that for N > 1, the vector 
(wfc^iv, 1 < A; < A^) follows the distribution Diriv(a) (a > 0). Then, the family 

N 

/i^ := ^Wfc,jv%,^ (37) 

k=l 

satisfies the LDP in A1i([xi)X2]) with speed (N) and good rate function is 
a}C{u\-). 

Proof: 

The proof uses heavily the classical representation ([7j) : 

{w,^,l<A;<iV}^i f ^ ^ ^ (38) 

where the F's are independent and 7a distributed. The modified measure 

N 



k=l 



has then independent weights and is easier to handle. We come back to the 
original measure with 

(TV) _ J^N 
/iTv(l) 

Now, recall that the cumulant generating function of the 7^ distribution is 
Lvir) 



alog(l — r) if r < 1, 
-00 otherwise. 
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and that its Cramer transform is (see for example [S]) 

. , fx — a — a log- if X > 0, 

l+oo otherwise. 

To understand the flavor of our proof, let us first assume that the sequence 
{Ck,N) is not random, and that the sequence (Cn) converges weakly to u. 
In that case, using the main Theorem of pO] we see that the sequence /I at 
satisfies the LDP with speed (A^) and good rate function : 

/(/^):=/ ijai9)-9)diy+ [ d/i (/i G A<([xi, X2])) , (39) 

where dfi = g du + {fi — g du) is the Lebesgue decomposition of /i with respect 
to z/. Further, as almost surely AtAr([xi, X2]) > and using the contraction 
principle (see [9] p. 126), we may conclude that {^n) satisfies the LDP with 
speed (N) and good rate function 

/(O = inf {/(/i) : /iGA^([xi,X2]), /i = /i([xi,X2])e} = mf7(&0- 

b>0 

Therefore, a direct computation yields 

1(0 = alC{u\0 . 

Let us now turn to the general case. We have to get rid of the randomness 
of the ^k,N- But, we show that, since the LDP with larger speed holds for its 
empirical measure, this randomness has no effect on the rate function. To 
get the LDP, we use a method inspired by the proof of the so-called Gartner- 
Ellis-Baldi theorem ([9] p. 157). Roughly speaking, it consists of two steps. 
First we compute the limiting normalized generating function and its convex 
conjugate function. Then we conclude by density of exposed points of this 
last function. 

To begin, let / be a continuous function on IR such that Lyof is bounded. 
Then, integrating first on the random variables Y we may write 



'AT 



Eexp/I^(/) = Eexp [ Ly (/(^fc.iv))] = Eexp (Ar£^(Lyo/)) . 



Now, from the assumption, the sequence of random variables {Cn{Ly° f)) 
converges to JlLyof) dv and satisfies the LDP with speed (A^''). Now fix 
e > and set 

A:={r^e -M+([xi,X2]) : |r/(Lyo/) - z/(Lyo/)| < e}, 
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and In = In, a + In,A'' with 

In,a := E [exp {NCn{Lyo f))l{c^^A}] • (40) 
Obviously, we have 

exp[iV (z/(Ly o/) - e)]¥ {Cn e A) < In,a < exp[N (z/(Ly o/) + e)] (41) 

and 

In,a^ < exp[iV||Ly o/ll^] P {Cn G A'=) (42) 
The LDP upper bound gives 

hm sup ^ log F{Cn eA')<- inf { (r/) ; r/ G A^} . (43) 

Since is good, its infimum on the closed set is reached and since it has 
a unique global minimum at u which is not in A^, the right side of (l43l) is 
negative. Since p > 1 we get limsup^ -i: log /iv,A'^ = —oo, and then 

limsup^log/iv < z/(Lyo/) + e. (44) 
To get the lower bound, we have 

liminf log In > liminf ^ log /at a > i^iLyo f)—e + liminf ^ logP (/Zat G A) 

N N N N ' - ^ -I > N N ^ 

From (113]) it is clear that limP (^at G A) = 1 so that 

liniinf ^ log /at > z/(Ly o /) - e (45) 

As ( 144]) and (145|) hold for every e > 0, we may conclude 

1 

lim — log/AT = z/(Lyo/) . 

The remaining of the proof is the same as the LDP proofs developed in [T7] 
(see also [ID]). ■ 
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5.2 Large deviations for spectral measures in the cir- 
cular and Jacobi ensembles 



We now state the LDP for the sequence of random probabihty measures 
ifJ'i^^) in the above frameworks. Notice that in the case /3 = 2, with the help 
of Theorem 14 ■![ the first result has yet been shown in |27] and the second one 
in [15] Theorem 2.3. In this last paper, the proofs use the so-called projective 
limit approach. We take here a completely different way. Indeed, we give a 
general proof by using a variation on the results developed in [HD] for random 
measures. As a matter of fact, we have shown in Subsection 15 . 1 1 1 he LDP for 
random measures with weights having a Dirichlet distribution. 

Theorem 5.3 

1. Assume that for any N, x Sn is endowed with the C/?Ejy distribu- 
tion. Recall that fii^^ is built as in ( f^) . Then, the sequence of random 
probability distributions {ni^^) satisfies the LDP with speed N and good 
rate function 

J(0 = |/C(At|0, (ee^i(T)) (46) 
where Xj is the uniform distribution on T. 

2. Let a,b > be fixed. The sequence {^^^^) under the J{P, a,b, N) 
Dirjv(/9/2) distribution satisfies the LDP with speed N and good rate 
function 

1(0 = ^/C(arcsine |0 (47) 

Proof: The proof of this theorem is a direct consequence of Proposition [521 
Indeed, from [20] it is known that the empirical distribution built on {6k) 
satisfies the LDP with speed N"^ whose rate function has a unique minimum 
in At. For the Jacobi ensemble, Hiai and Petz ([22] Section 2) proved the 
LDP with speed A^^ for the empirical measure built on (x^). Moreover, they 
have shown that - with our assumptions - the unique minimizer of the rate 
function for this LDP is the arcsine distribution. 

Remark 5.4 The same conclusion holds true under the probability measures 
of TheoremU. 
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Remark 5.5 The assumptions about the supports are crucial. In a compan- 
ion paper IT^ we study important models in which they are not satisfied. For 
the GUEtv, the support of Cn is not included in a fixed compact set and the 
limiting distribution is the semicircle. In this example, we find an additional 
term in the rate function. For the J {(3 , tiN , T2N ; N) , the support of Cn is 
included in [0, 1] but the limiting distribution is supported by a compact subin- 
terval of (0, 1). It is called the Kesten-MacKay distribution and generalizes 
the arcsine distribution. In ITO^ . we prove the LDP but the explicit form 
of the rate function is only conjectured. The same is true for the Laguerre 
ensemble. 

Remark 5.6 Other discrete random measures with Dirichlet weights appear 
in the literature. Let a be a positive measure on a compact space K . A 
Dirichlet process fi of parameter a is a probability distribution denoted by 
©(a) on Aii{K). It is such that for any measurable finite partition of K, 

{Ai, ■ ■ ■ ,An), we have (/i(v4i), ■ ■ ■ , fi{AN)) = DirAr(a(y4i), ■ ■ ■ ,a{AN)). 

If a is the finite discrete measure 9 J2k=i ^Vk where yi, - ■ ■ ,yN & K, we have 



where (Yi, ■ ■ ■ , Y/v) = DiTjsf{9) . This means that the random measures studied 
in the previous sections are mixed Dirichlet processes. The mixing distribu- 
tion is the law of the random eigenvalues. Ganesh and O'Connell (fTEp study 
V{a + Xlili ^Vi) when J2iLi converges to some measure v and when 
a is a measure whose support is K . They show that /x satisfies the LDP with 
speed N , and good rate function /(■) = /C(z/|-). In ^ there is an extension 
with an infinite number of random locations, (see also IS^ for the connection 
with the Poisson- Dirichlet distribution). 

5.3 Relation with spherical integrals 

In the unitary model, the direct computation of the hmiting cumulant 
generating functional of leads to a spherical integral. Indeed, if G C (T) 
(the set of all continuous function on T), we have 




(48) 



(^)(^) = (ei,^(f/)ei) = tr(G^yDjf)\/*) 
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where D^^^ := diag(v9(e*^^), ■ ■ ■ , ^^(e*^^)) , (the e*^* are the eigenvalues of U), 
V G U(A^) and := diag(l, 0, ■ ■ ■ ,0). Consequently, the Laplace transform 
of yU„ is, for Lf e C{T), 

Eexp (iV/.f)(^)) = E,/(J)(G^,D(f)) (49) 

where 

iS\G^,dI^^) =( [ exp[iVtr(G;vV^I^JfV*)] dV 

\JV(N) 

(Eg denotes here expectation with respect to the variables 9i, . . . , Ojy). It is 
interesting to notice that this last spherical integral over the unitary group 
can be expressed as a hypergeometric function with two matrix arguments 
(see for example [31] p. 97). 

In [19] and [6] it is proved that N^^ log I^^\Gn, Dyy^^) has a limit as 
N ^ oo, as soon as the empirical spectral distribution of d!^^ converges 
weakly. More precisely in [19] Th.6 (see also [6] Th. 4) it is proved that 

limilog/iJ)(G^,D(f)) = F^(l), (50) 

where -F(^(l) is defined in the following way. Let be the limit of the 
empirical spectral distribution of (in our case it is the image by (p of 

At). Further let (pmin '■= min^gf V^(2;) and v^max := max^gx ¥^(2;). The Stieltjes 
transform of z/^^ is defined for x G {—00, v?min) U (v^max, +00) by: 

We have that := lim^i^^^^ H^{x) < < := lim^|<^^^^.^^ if^(x) and 
the range of is {H^,0) U (0,i/^). ([19j Property 9). The Voiculescu 
i?-transform is the function satisfying, for all y in the range of H^, 

H^(^R^{y) + ^^ = y. (52) 

As a result ([TU] Theorem 6), 

F^(l) = v{l)- [\ogil + v{l)-y,iz)) 

V{1) = <<^max-l ifl>i^i (53) 
[V^min-l ifl<ift 
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Now, the LDP for the empirical distribution built on 6 



1) 



6n holds with 



speed N^. So that, A^-^ InEexp {N /j,];^ > {cp)) and A^-^ log^ (Gjv, ^if ) have 
the same limit (to show this claim just use the continuity of spherical integrals 
proved in [28] Prop. 2.1 

recovered by taking the supremum in y9 G C(T) of fi{(f) 
H e Mi{T). Setting g{e'^) = (/?(e^^) - v{l), we have 



The rate function of the LDP for /ii '' can be 



(N) 
w 

where 



fi{^) - F^{1) = / g{z)fi{dz) + / log(l - g{z)) Xjidz) . (54) 
Taking the supremum in g E C(T), we recover the well-known duality formula 



sup 

96C(T) L^T 



g{z)d^{z) + / log(l - g{z)) Xji^z) 



/C(AT|/i) 



References and several consequences of this formula for the associated mo- 
ment problem may be found in [15]. 
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